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Abstract 



^ (— I The most general combination of couplings of fermions with external potentials in 1+1 

O |' dimensions, must include vector, scalar and pseudoscalar potentials. We consider such a mix- 

ing of potentials in a PT-symmetric time-independent Dirac equation. The Dirac equation is 
mapped into an effective PT-symmetric Schrodinger equation. Despite the non-hermiticity 
of the effective potential, we find real energies for the fermion. 
PACS numbers: 03.65. Ge, 03.65.Pm, 03.65.Fd 



1 Introduction 

> 
O 

^r) [ In the last ten years the, so-called, PT symmetric systems introduced in the seminal paper by 

0^ ■ Bender and Boettcher [2] have attracted very much attention. In fact, there are many works 

devoted to develop and understand this new kind of situation pQ 119) . They consist of non- 
Hermitian Hamiltonians with real eigenvalues, which however exhibit parity and time-reversal 
symmetries. The one-dimensional time-independent Schrodinger equation is invariant under 
space-time inversion. In addition, there exist other classes of Hamiltonians with real spectra 
without being PT symmetric, as can be seen, for instance, in the references [HIE]; and systems 
where the PT symmetry is spontaneously broken with complex energy eigenvalues [7j. The 
problem of non-Hermitian time-dependent interactions which does not exhibit PT symmetry 



but still admits real energies is considered in the references [9]-|12|. 

More recently, some problems of relativistic fermions interacting with non-Hermitian poten- 
tials of scalar and vector natures have been reported in the literature [T].[19j-[29]. In general, 
it has been shown that for some configurations of those non-Hermitian potentials, the Dirac 
equation admits real energies. In 1+1 dimensions the Dirac problem is easily mapped into 
a Sturm-Liouville problem or, in other words, into a time-independent Schodinger equation 
with real or complex potentials whose bound-state solutions present real energy eigenvalues. 
Some authors [l3]-[18] have investigated the Klein-Gordon equation and Dirac equation in the 
context of PT symmetry and pseudo-Hermiticity. Mustafa [13] studied the exact energies for 
Klein-Gordon particle and Dirac particle in the generalized complex Coulomb potential. Zno- 
jil [2] analyzed the Klein-Gordon equation presenting a pseudo-Hermiticity. Egrifes et al |15j 
investigated the bound states of the Klein-Gordon and Dirac equations for the one-dimensional 
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generalized Hulthen potential within the framework of PT-symmetric quantum mechanics. In 
|17j . the authors have investigated the bound states for the PT-symmetric versions of the Rosen- 
Morse and Scarf II potentials in the Klein-Gordon equation with equally mixed vector and scalar 
potentials. Sinha and Roy [18] investigated the one-dimensional solvable Dirac equation with 
non-Hermitian scalar and pseudoscalar couplings, possessing real energy spectra. 

The coupling with scalar potentials in the Klein-Gordon and Dirac equations can be seen as 
a position-dependent effective mass. In the relativistic ambiance, the ordering ambiguity of the 
mass and momentum operators, which is present in the non-relativistic one, should disappear. 
Nevertheless, there are difficulties to define consistently fermions and bosons, whenever one takes 
into account space-time dependent masses. This happens due to the fact that physical particles 
in quantum field theory must belong to an irreducible representation of the Poincare algebra 
[301 13 lj . One should be able to find generators specifying the particle properties, usually its 
mass and helicity. However, it is quite hard to accomplish with this task in the case of spatially 
dependent masses. Thus, one should keep in mind that all of these usually thought as relativistic 
equations for position-dependent masses should be taken as effective equations. In this regard, 
Alhaidari [32] studied the exact solution of the three-dimensional Dirac equation for a charged 
particle with spherically symmetric singular mass distribution in the Coulomb field. Vakarchuk 
|33j investigated the exact solution of the Dirac equation for a particle whose potential energy 
and mass are inversely proportional to the distance from the force center. In [26] . the authors 
considered the smooth step mass distribution and solved approximately the one-dimensional 
Dirac equation with the spatially dependent mass for the generalized Hulthen potential. In |34j . 
the authors investigated the exact solution of the one-dimensional Klein-Gordon equation with 
the spatially dependent mass for the inversely linear potential. Therefore, it is of considerable 
interest to investigate the solution of the effective mass Klein-Gordon and Dirac equations with 
non-Hermitian complex potentials with real energy spectra. 

One interesting problem that has been tackled in this context is the Dirac equation in 1+1 
dimensions in the presence of a convenient complex vector potential plus a real scalar potential 
[I], wherein the scalar potential plays the role of a position dependent mass. Here, we will show 
that one can realize a more general system of massless fermions in two dimensions interacting 
with a mixing of complex vector, scalar and pseudoscalar potentials. Although complex, the 
scalar and pseudoscalar potentials are responsible to open a mass gap for the fermions. Such a 
scenario might be important for some condensed matter systems, where the electrical conduction 
is essentially one-dimensional. The scalar and pseudoscalar potentials can be thought as defects 
in the lattice and the electrons are also subject to a background potential of vector nature due 
to the ions in the lattice. 

Our purpose in this work is to show that the Dirac equation in a two-dimensional world can 
still have real discrete energy spectrum and supports fermion bound-states when a convenient 
mixing of complex vector, scalar and pseudoscalar potentials is considered. We call attention 
to the transformation of the potentials under parity in order to have the PT-symmetry in the 
Dirac equation. Particular configurations of those potentials are worked out in some detail. The 
approach here is the mapping of the PT-symmetric Dirac problem into a naturally PT-symmetric 
effective Schrodinger equation. 
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2 The time-independent Dirac equation in 1+1 dimensions 



We consider here the 1+1 dimensional time-independent Dirac equation for a massless fermion 
under the action of a general potential V. It is written as 

HV(x) = EV[x) , (1) 

H = cap + V, (2) 

where E is the energy of the fermion, c is the velocity of light and p is the momentum operator. 
a and (3 are Hermitian square matrices satisfying the relations a 2 = (3 2 = 1, = 0. 

The positive definite function |^| 2 = ty'ty, satisfying a continuity equation, is interpreted as 
a position probability density. This interpretation is completely satisfactory for single-particle 
states. 

We set V to be 

V = (5M{x) +f3-f 5 P(x) + V{x) +(5aA{x) (3) 

where M(x) is a scalar potential, P{x) a pseudoscalar potential and V{x) is the time-component 
of a Lorentzian 2-vector potential, whose space component is A{x). The space component of the 
2-vector potential can be eliminated by a gauge transformation without affecting the physics. 
Once we have only four linearly independent 2x2 matrices, the structure of coupling in V is the 
most general one can consider in the time-independent Dirac equation in one space dimension. 
In terms of the potentials the Hamiltonian ([2]) becomes 

H = cap + V(x)+f3M(x)+(3-f 5 P(x), (4) 

where 7 5 = — ia. An explicit expression for the a and j3 matrices can be chosen from the Pauli 
matrices that satisfy the same algebra. We use f3 = a\, a = cr 3 , and thus (3j 5 = —02- The 
equation (pQ) can be decomposed in two coupled first-order differential equations, for the upper, 
and lower, components of the spinor In a simplified notation and by using 

the natural system of units h = c = 1, we have 

+ ifj/_ + MV>+ + Vip- - iPip+ = Etp- . (5) 

where the prime stands for the derivative with respect to x. 

In the case that all the potentials are real functions, the Dirac equation is Hermitian and 
invariant under space-reversal (parity) transformation. We recall that the parity transformation 
is an improper Lorentz transformation and that the spinor in one frame is constructed from 
the spinor in the other frame by means of the relation ^(x,t) = S^(x,t) = e lS /3^>(x, t), with 
x = —x and 5 an overall constant phase factor. Moreover, under parity transformation, M{x) 
and V(x) do not change and P(x) changes its sign. The matrices must transform as 5'~ 1 /35 = (3 
and S~ 1 aS = a. 

The same transformations could be used even in the case that the potentials are complex, but 
if we want the Dirac equation invariant under the combination of parity and time-reversal trans- 
formations this issue becomes trickier with a non-Hermitian Hamiltonian. This is because the 
time-reversal transformation implies that T(i)T _1 = — i and that the potentials in the Hamilto- 
nian are complex. Thus, althougt the time-reversal does not change each part of the potentials, 
since they are time-independent, it changes the relative sign between the real and the imaginary 
parts of the potentials; as a consequence, the imaginary part of each one of the potentials must 
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change under parity in the reversed form of its real part, in order to have the Dirac equation in- 
variant under the combination of parity and time-reversal transformations. In summary, in order 
to have PT-symmetry even when the potentials are complex, the imaginary part of the vector 
and scalar potentials must change their signs under parity, whereas the imaginary part of the 
pseudoscalar potential does not change. The spinor in the time-reversed system is obtained from 
the spinor ^(x,t) by means of the following the transformation ^q-(x, t) = T^(x, t) = T^/*(x, t), 
with t = —t and T a square matrix such that T~ 1 (3*T = (3 and T~ 1 a*T = —a. Then T must 
commute with (3 and anti-commute with a, that is T = (3. 

All the examples of potentials we are going to deal with follow the transformation rule given 
above. 



3 The effective PT symmetric problem 



Whenever one considers only the coupling either with the scalar or the pseudoscalar potential, 
the differential equations can be uncoupled in such a way that both components of the spinor 
satisfy second-order differential equations, similar to each other and to the Schrodinger equation. 
By including the coupling with the vector potential this is no longer possible. Although, it is 
possible to show that one of the components obeys a kind of Schrodinger equation, and the other 
component is given in terms of the previous one, as was done in the references in pQ. From now 
on, we are going to follow that approach. By applying the space derivative in the first of the 
equations ([5]) we have 

- + + A' + 4>_ = B^' + + (6) 

where we have defined A± = A±{x) = M(x)±iP(x) and B = B(x) = E — V(x). By substituting 
in the above equation the expressions for ip'_ and ip- taken from equations ©, we obtain the 
following equation for the upper component 



BAj 



B' + i(A + A^ - B 2 



0. 



and the equation obeyed by the lower component can be rewritten as 



1 



— (#V + Bip^ 



We notice that by means of the redefinition of the upper component 

the first-derivative term can be eliminated and x satisfies the differential equation 



- IX" + 



BA\ 

~a~ 



B' + i (A+A- - B 2 ) + 



2 Al 



x = o, 



which, in its turn, can be written as a Schrodinger, equation 

-x" + V effX = E 2 X , 

where 



V eff =E[2V + i 



A' 

'A+ 



iV' + i 



A+ 



V + 



A+A- — V + 



A+ 



A+ 



(7) 



(8) 



0) 



(10) 



(11) 



(12) 
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Equations (f7|)- ([T2l) constitute the essential set of equations we are going to deal with through- 
out this work. Our intention is to find convenient configurations of the potentials for which 
the equation (llip presents normalizable solutions and positive eigenvalues or, in other words, 
bound-state solutions corresponding to a discrete energy spectrum, with E 2 > 0. Naturally, the 
solutions for the lower component that satisfies equation ([8]) must also be normalizable. 

In the literature we find several studies of the Dirac equation in a two-dimensional world 
with convenient mixing of scalar plus vector potentials and pseudoscalar plus vector potentials. 
As far as we know, most of them assume an intrinsic relation between the scalar and vector 
potential, or between the pseudoscalar and the vector potential. A typical case is the one which 
deals with a constraint of the type M(x) = nV(x), where \k\ > and M(x) is a binding 
potential [35j. There are some cases in four dimensions for which such convenient relation leads 
to analytical and exact solutions for the Dirac equation. Here, we are not going to escape from 
a constraint among the potentials. On the contrary, we follow the same procedure carried out 
in the references [I] by constraining the vector potential to obey the more general relation 

v = i(^±) + v = i( M ' +ip, ) + y. (is) 

2 \A + J 2\M + iP ) { J 

With this relation the Schrodinger-like equation (jlip is given by 



X + 



M + P —V + 2EV -iV 



X = E 2 X . (14) 



The effective potential written explicitly in terms of the real and imaginary parts of the 
potentials 

V ef f (x) = [m 2 - M 2 + Pi - P 2 - V r 2 + Vi + 2EV r + Vty + 

+2i (M r Mi + P r p - V r Vi + 2EVi - vf) (15) 



is complex, but the effective Schrodinger equation exhibits PT symmetry if M r , V r and Pi 
remains invariant, whereas Mj, Vi and P r changes their signs under parity transformation, as 
was stated in the second section. Naturally we are thinking of the cases for which the energy in 
the Dirac problem is real (positive or negative). 

4 Real spectrum with non-Hermit ian potentials 

As our first example, we choose a set of potentials which have a quadratic dependence on the 
spatial variable and, in order to deal with a complete set of exact eigenstates, we impose some 
restrictions over the potential parameters. In this first example we choose 

M{x) = a x 2 , P{x) = i(a lX 2 + h), V(x) = a 2 x 2 (16) 

In this situation, the vector potential is given by 

V = 1 [ aox>-(a lX * + bl) \ + ^ ■ (17) 
Consequently the effective potential is 

Veff = (oq -a\- al)x 4 + (2Ea 2 - 2aih)x 2 -b\- 2ia 2 x. (18) 
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As we asserted, we want to deal with an exactly solvable model. We restrict the potential 
parameters in order to eliminate the quartic term, obtaining the constraint ao = yaf + af. 
This lead us to the following effective PT-symmetric driven harmonic potential 



Veff = Q 2 x 2 - 2ia 2 x - b\ 



(19) 



where Q, 2 = 2Ea 2 — 2a\b\. By means of the transformation x = y + 2ia 2 /fl 2 we obtain the 
following differential equation for x(y) 



E 2 + b\ 



X, 



(20) 



Q(E) 2 

that is a quantum harmonic oscillator of unity mass. Consequently the energy E of the rela- 
tivists system is obtatined from the following equation 



E 



-b\ + 



n(Ey 



+ (2n + 1)0 (E) . 



The above equation is considerably simplified if we set b\ = 0, that is 

E 2 = ^ + (2n + l)^2E~V 2 , 

which can be written as 



(£ 3 -|)=(2n + l) E^/2E~V 2 , 



and, once squared, lead us to 



e 2 - a 2 



1 + 2(2 n + iy 



e + 



0, 



where we have defined e = E 3 , leading finally to the solutions 
E± 



y \l + 2 (2 n + l) 2 ± \l 1 + 2 (2 n + 1) 



(21) 



One can note that it is not necessary the presence of V{x) in order to have an effective 
PT-symmetric Schrodinger equation with a complex effective potential. As a matter of fact one 
can see that the dependence of V on the coordinate implies in an effective potential explicitly 
dependent on the energy of the Dirac problem and, as a consequence one will reach invariably 
an intricate transcendental equation for the energy. Henceforth, we set V a real constant. 

As an example, we suggest the following configurations for the scalar and the pseudoscalar 
potentials 

M(x) = ioj\x + mi, P(x) = uj 2 x + im 2 (22) 

with u>i, mi, uj 2 and m 2 real constants. The vector potential, with V = 0, is complex and given 
by 

v( ) = (m ~ ^2) [(mi - m 2 ) + i(m + uj 2 )x] 
{ ' " 2[( Wl +w 2 )¥ + (mi-m 2 ) 2 ] ' 1 j 

and the effective potential in the Schrodinger equation obeyed by x( x ) 

V e ff(x) = uj 2 x 2 + 2i(u\m\ + uj 2 m 2 )x + (m 2 — m 2 .), (24) 
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with ui 2 = lo\ — uj\ , is PT-symmetric. 

By defining the variable y = x — i(wimi + W2TO2)/w 2 we obtain 



dy 



^ + uj 2 (y 2 + A 2 ) X = E 2 X , (25) 



whose non-normalized eigenfunctions and the eigenvalues are respectively given by 

Xn (y) = e-"y 2 l 2 H n (^y) (26) 

and 

E 2 =w(2n + 1) + AV, (27) 

where H n {^/uJy) is the Hermite polynomial of degree n and A 2 = ( CJim ' 2 Jii a;2mi ) 2 - Consequently, 
the upper component of the Dirac spinor can be written as 

V>+,„(y) = N^(y 2 + \ 2 ) l l A e-^y 2 -^H n {^y), (28) 

where 6{y) = tan^ 1 (y/A), and N is a normalization constant. By recalling the equation ([8]) and 
the expression for the vector potential, V(x) = iA' + /2A + , we get the following equation for the 
lower component of the spinor 

*- w -7m(^ +Ex )- (29) 

which can be written in terms of Hermite polynomials as 

1>-, n (y) = ^ {y 2 N + X 2 ) i/4 e ' hiUJy2+l6{y)) [(± 1^1 - ^ H n{^y) + i2n^H n ^(^y)] . (30) 

The upper (lower) sign of the factor ± \E n \ in the expression above stands for the positive 
(negative) energy solution and \E n \ = y / w(2n + 1) + \ 2 u> 2 . The normalization constant N is 
obtained by the definition of the norm of the eigenstates 

hoo r+oo 

|2j„ _ / (U, A..\|2 , I,,/, f„,M\2, 



\* n \ 2 dx = / (\^ n (y)\ 2 + \^ n (y)\) 2 dy = l, (31) 

oo J — oo 

which can be written in terms of the variable u = \fjjy as 
/■+oo e - " 2 

2 l^l 2 / du rzTTa M/2 { + A ' w ) + ( 2n + X )] " ^nH n ^(u)H n+1 (u)} = 1, 

(32) 

where we have used the recursion relation 2uH n (u) = 2nH n ^i(u) + H n+ i(u). It is very hard 
to check the convergence of the above integral, for all values of n, even though each of them 
converges, since the weight function decays faster, for x — > ±oo, than any polynomial appearing 
in the numerator of the integrand. We furnish below the normalization constant for the two 
lowest energy levels, 

A?o = {e" A2w/2 (l + 2A 2 w)i^o(AV2) + v / ^C/(l/2,0,A 2 u;)} 
Ni = {e^/ 2 [2(1 + \W)K {\ 2 lj/2) + 2A 2 w(l - \ 2 u)K 1 {\ 2 u J /2)] + 
+V^r(l + 2A 2 w)?7 (1/2, 0, A 2 tj)}~ 1/2 , 
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where Kq(z/2) and K\(z/2) are modified Bessel functions of the second kind and U(a,b,z) is 
the confluent hypergeometric function. 

As a third example we take the scalar and pseudoscalar potentials given respectively by 

M = iM x tanh(/xx) + M , P = Pi tanh(^x) + iP , (33) 

with Mi, Pi, Mo, Pq reals and constants. Thus, the vector potential is also complex and, with 
V = V , is given by, 



fx (Pi + Mi) [(Mp - P ) - i{Mx + PQ tanh(/ix)] sech(^r) 
2 (M -P ) 2 + (Mi + P 1 )2tanh 2 (^) " l " 



and the effective Schrodinger equation (|14|) can be written as 

- X" - = £X , (35) 

where 

V e ff(x) = -{Pi - M^sech 2 ^) + 2i(M M 1 + P Pi) tanh(^x) (36) 

and 

e = (E - V ) 2 - (A 2 - M 2 ) + (P 2 - M 2 ) . (37) 

This effective potential is a generalization of the real Rosen-Morse II (hyperbolic) |38|, [39] and 
was already considered in the literature concerning PT symmetry in non-Hermitian Schrodinger 
problems [5] , [6] . The eigenvalues of the effective energy can be conveniently written as 

e n = -n 2 (a 2 n -b 2 n ), (38) 

where 

a n = (1/2) 



1 + (2?7i//i) 2 - (2n + 1) 



and b n = a n , (39) 



with Uf = P{-M{, n = M M 1 +P P 1 and n = n = 0, 1, 2, ... < s = (1/2) 
Thus, the energy eigenvalues of the Dirac problem are given by the following relation 



1 + {2U 1 /fi) 2 - 1 



E n = V ± ^Uf - C/ 2 - ^{al - 62), (40) 

where Uq = Pq — Mq. In order to have a real spectra, the relation XJ\ — Uq > A i2 (a 2 — fr 2 ) 
among the parameters of the potential must be imposed. It worths mentioning that this kind of 
restriction is not a consequence of a complex potential in the effective Schrodinger equation. It 
appears, instead, due to the kind of coupling and to the configuration of the potentials we have 
considered in the relativistic problem. 

If the potential supports n max + 1 bound states, then we can say that s = n max + 5, with 
< 5 < 1. We have observed that, for Q 2 = /i 4 <5 4 , the factor /U 2 (a 2 — 6 2 ) decreases in the 
range [0,ra max ]. The maximum value achieved by ^ 2 (a 2 — 6 2 ) is fi 2 {s 2 — 5 4 /s 2 ) . Therefore, the 
relations 

5 4 



J 2 = fj, , Ui — Uq > fl 



("max + S)' 



(n max + 6) 



(41) 



are necessary to have a real spectrum. In the case we have only one bound state for example, 
corresponding to n max = 0, we obtain M Mi + P Pi = [i 2 5 2 , P 2 - M5 2 < P x 2 - M 2 = fi 2 5(5 + 1). 
In general the upper component of the spinor is given by 
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il>+, n (z) = N+\A+(z)\ 1/2 e ia ^ 2 {l - z 2 ) a " /2 



1-z 
l + z 



p(a n +b n ,a n —b„) 



(42) 



where, z = tanh(/irr), P^ 1 ^' \z) are the Jacobi polynomials of degree n |42j . On its turn the 
lower component can be written as 



V>-n(*) = N_\A + (z)\-^ 2 e~ ia( - z ^ 2 fi(l- z 2 ) a ^ 2 



1 - z 
l + z 



d 



(E n - a n z - 2b n - v )p( an+bn > an - bn) + (1 - z 2 )_pK+^-M 



(43) 



In the expressions (02} and (03]) we have |^+(z)| = y/{M - P\z) 2 + {M\z + P ) 2 and a(z) 

tan- 1 ((M 1 z + P )/(A/ -i :, i^). 

The normalization constants are found by the definition of the norm, namely 



+1 



dz (i-z 2 )- i (i^+,nr + iV'-,nr) = i. 
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5 Final considerations 

We have extended the problem of fermions in 1+1 dimensions interacting with complex potentials 
by considering the most general coupling of fermions which is possible in the time-independent 
Dirac equation in two dimensions. The potentials are considered to be of scalar, pseudoscalar 
and vector natures. We have shown that one can map the Dirac equation with all those complex 
potentials into a PT-symmetric effective Schrodinger equation, which can be exactly solved and 
present real spectra. Indeed, the PT-symmetry is also observed already in the Dirac equation 
when the imaginary part of each one of the potentials transforms, under parity transformation, 
in the reversed manner that the real part of the potential. In this way, the PT-symmetry seems 
to play a fundamental role for the reality of the spectrum. There still remains the open question 
about pseudo-hermiticity of the Dirac operator with the mixing of potentials we have considered. 

Another remarkable point of those examples, besides being exactly solvable, is that they bind 
fermions despite the presence of the vector potential. Inhomogeneous vector potentials give rise 
to electric fields which are responsible for pair production. Based on what happens when we 
have a convenient mixing of scalar and vector potential [35], we think that the presence of the 
scalar and pseudoscalar potentials, although complex, are responsible for giving an effective mass 
to the fermion and, as a consequence, the threshold for pair production becomes higher than 
the energy that the electric field can supplies to the fermion. Confinement of fermions with a 
vector potential, as in the second example treated here, is a peculiar case. Indeed, one can check 
that the relativistic spectra we have found is very similar to that of the "Dirac oscillator" in 
1+1 dimensions [33]. In that reference a massive fermion is coupled to a linear potential. That 
coupling, being of a pseudoscalar type, is the two-dimensional version of the four-dimensional 
coupling of an anomalous magnetic moment with an electric field which emerges in the "Dirac 
oscillator" [33J-[37], and that can be seen as a mechanism for confinement of fermions in 3+1 
dimensions. 

We have noted that in the absence of the scalar and pseudoscalar couplings, the effective 
equation for the upper component of the spinor is given by 

-t/j'l +[2iV'-(E- V) 2 )] = 0, 
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that is, we have a PT-symmetric effective Schrodinger equation with a vector potential which 
is real and even under parity transformation. This is a situation where the PT-symmetry 
becomes quite important, because it opens the possibility of having real spectra in a model 
which, otherwise, would be discarded due to the fact that the effective Schrodinger equation 
corresponds to a non-Hermitian operator. This opens the promising perspective of finding 
exactly solvable real potentials where this feature could be accomplished. 
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